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BALLOON TEST PROJECT 
COSMIC RAY ANTIMATTER CALORIMETER (CRAC) 

J.C. Christy, G. Jhenain, P. Goret, J. Jorraiid, 
P. Masse, P. Mestrea.u, N. p etro”., A. D o^ir. 
Astrophysics Service (CEN-Saclay) 


I - INTRODUCTION 


/l* 


The presence of antimatter in galactic cosmic radiation was 
recently established (1), (2), (3). The table below summarizes 

the characteristics of these three experiments. 


Ref . 

Energy 

Geometr . 

t of p 

Type of 

P/P 


(Gev/n) 

Factor 
(cm2 . sr ) 

de tec ted 

detector 

measured 

(1) 

3-6 

200 

30 

Magnet 

5xl0~ 4 

(2) 

2-5 

1 

2 

Magnet 

6xl0 _4 

(3) 

< . 3 

120 

14 

Calori- 

meter 

-4 

2x10 


Figure 1 shows the p/p ratios measured as a function of 
energy. The curves correspond to the predictions of the various 
cosmic radiation propagation models in the galaxy, on the hypo- 
thesis that there is no primary antimatter component. We see that 
the two high-energy measurements are fairly close to the predic- 
tions, whereas the low energy point is on several orders of magni- 
tude Above the expected value. This result could indicate the 
presence of a primary antiproton component. Yet, this conclusion 
remains very hypothetical for two reasons: 

i- the detection method used for low energy conditions renders 
the result very uncertain; 

ii- not a single antihelium nucleus has been detected yet 
(He/He < 2x10 5 according to ref. (3)). 

* Number s in the margin indicate pagination in the original text. 
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Another explanation recently suggested concerns the n-n 
oscillation predicted by the theories of high unification (4). 

The purpose of the CRAC experiment is to measure the flux of anti- /2 
matter (p and He) in the 200-600 Mev/n energy range-, with the 
following characteristics: 

-a high geometric factor (420 cm2.sr) 

-a high detection sensitivity 

-the presence of contamination due to spurious events <10 6 

These three factors combined should make it possible to de- 
tect antimatter at the p/p=10 4 level. 

No other experiment on antimatter in cosmic radiation is 
being considered by other laboratories. 

The second objective of the CRAC experiment is to measure the 
isotopes of light elements (Z<8) between 600 and 1000 Mev/n. The 
measurements performed until present pertained to low energy (<150 
Mev/n) and their interpretation is obscured by i) solar modulation 
which is very substantial at low energy and ii) our inadequate 
knowledge of the effective sections which also vary rapidly with 
energy. The measurement proposed here overcomes these shortcom- 
ings for the most part. The scientific benefit of these isotopic 
measurements is that they make it possible to determine the RC 
trajectory distribution in the galaxy corresponding to the long- 
est trajectories. 

II - p Detection Principle 

When an antimatter nucleus interacts with ordinary matter, 
this causes annihilation associated with the emission of second- 
ary pions. Figure 2 shows a diagram of a p-p annihilation. On 

the average, one annihilation produces 3 charcre Dions (tt + , tt ) and 
two neutral pions (tt°K The tt 0, s promptly disintegrate into 2y 
which may produce a pair (e + ,e-) , or a Compton electron, due to 
its interaction with surrounding matter. The energy balance of 
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The energy balance of the annihilation is: 

Available energy 1.9 Gev + kinetic energy 

Pion mass 0.7 Gev 

Total kinetic energy 
available for the 

pions 1.2 Gev + kinetic energy 

We thus see that the mean kinetic energy of a pion is about 
240 Mev. Figure 3 shows in more detail the energy distribution 
of the annihilation pions. We thus see that the annihilation pro 
ducts (TT,e+e-) are ultrarela tivistic particles. We will there- 
fore use a Cerenkov detector to detect the annihilation products. 
In our case, there is another advantage to using a Cerenkov de- 
tector:.; it is insensitive to slow moving particles. Moreover, 
as slow.; moving products are abundant at the time of the nuclear 
interactions within the Cerenkov itself, the background noise is 
thus considerably reduced. 

Ill - The Telescope 

The telescope proposed is shown in figure 4. It is made up 
of 2 modules: 

-module 1 has 5 functions: 

1) to trigger a coincidence signal between Tl and T2 

2) to identify the particle charge by means of 6 scintil = 
la tors Tl, T2, SI, S2, S3, S4 

3) to measure the flight time between Tl and T2 so as to 
identify the propagation direction 

4) to measure the velocity of the incident particle in the 
Cerenkov Cl 

5) to measure the trajectory using a hodoscope Sl,S2,S3,S4. 

-module 2 is a calorimeter made up of 5 Cerenkov C2 detector 
(1 to 5) . It is designed' to detect the Cerenkov light emitted by 
the annihilation products. 
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An anticoincidence protecting module 1 is provided to quickly 
discard events beyond coincidence . 

a- Flight time: 

This is the critical point of the experiment. There actually 
exists a large flux of rising particles either coming from the 
atmospheric albedo, or from nuclear interactions in the calori- 
meter. These particles are not only interesting, but are also 
capable of simulating an antiproton signal, if necessary. It is 
therefore necessary to discard these particles with 100% effi- 
ciency. In our case, the flight time between the falling parti- 
cles and the rising particles is < 6 ns. There are ..actually flight 
time devices having resolutions on the order of o=0.2 ns (ref. (5) 
for example) which would give us a comfortable separation of more 
than 3 0 sigmas. 

b- The Hodoscope: /4 

There are 4 layers (2X and 2Y) of scintillator rods (50x5x3 
cm) seen by a PM at each end. The angular resolution is < 2.0 
degrees. The advantage of such an hodoscope is its 100% proton ef- 
ficiency and its relative simplicity. The spatial resolution is 
amply adequate for correcting the Cerenkov signals of the chart 
effects which is the main cause of fluctuations. To detect anti- 
matter, the angular resolution is of little significance, the 
hodoscope being used in this case (excluding determination of the 
charge) only to contain the events in the calorimeter geometry. 

c- Identification of the Charge: 

It is well known that scintillators have an excellent reso- 
lution for low charge particles which we are concerned with here. 
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d - Cerenkov Radiators: 


The various Cerenkov radiators used in the experiment are 
all designed according to the same principle. They are made up 
of a diffusing box filled with a mixture of purified water & shift- 
er (2-amino-6 acid di-potassium salt, disiilfonic 3 -nap thtalene 
abbreviated amino-G salt) . The boxes are round with dimensions 
<j)=8 0 cm - h=20 cm for the 5 Cerenkov of the calorimeter and (J> = 60cm 
- h=15 cm for the Cl Cerenkov. The water in each box is seen dir- 
ectly by 12 PM, the PM window being in contact with the water. We 
thus gain a factor on the order of 3 to 4 on the number of photo- 
electrons with respect to a PM viewing the water through a trans- 
parent window. Numerous Cerenkov with wa ter+ shif ter were used in 
in the pass (6), (7). The addition of shifter has a dual effect on 
the primary Cerenkov light: i) the UV radiation which is more sig- 
nificant is absorbed by the shifter and reemitted at 400 run, i.e. 
in the photocathode efficiency peak and ii) the re-emission is 
isotropic which effaces the directivity of the primary emission. 

We thus gain a factor of about 4 on the number of photoelectrons 
by adding the shifter. Furthermore, an interesting property of 
water is the absence of scintillation, as shown in figure 5. The 
water ahs a refraction index of n=1.33 which gives an energy thresh 
old of 480 Mev/n. 

Based on the existing data on the Cerenkov with water, we 
may estimate the number of photoelectrons which may reach 
radiators such as those proposed here. The parameter which inter- 
ests us is the number of photoelectrons collected per cm of water 
for a charge particle Z=1 and velocity v/c=l. This parameter is 
shown in figure 6 as it is measured in the reference (6) as a 
function of the ratio of photocathode surface to the total internal 
surface of the radiator s-S (PM) / (total) S. We have also shown in 
this figure the number s corresponding to our radiators. We see 
that the Cl Cerenkov could collect 38 p . e/cm and the calorimeter 
Cerenkovs 25 p.e/cm. 
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Utilization of the Cerenkov water radiators thus appears 
to be quite appropriate for our present needs. One particle 
which passes the T1-T2 coincidence will trigger Cl or otherwise, 
depending on whether its kinetic energy per nucleon is less or 
not than the threshold energy of 480 Mev/n. The particle then 
penetrates into the calorimeter while decelerating. If the par- 
ticle is a proton it will give a calorimeter signal which is de- 
finitely lower than the Cl signal. Conversely, an antiparticle 
will become annihilated in the calorimeter by producing a signal 
which is much higher than the Cl signal. We will note that the 
antiparticles may become annihilated in flight or at rest after 
the deceleration. This fact affects the characteristics of the 
antimatter signal only very slightly. 

IV - Simulation of the Experiment 

A simulation of the annihilation process was performed by 
Monte-Carlo for the purpose of studying the telescope performances 
The simulation program handles: 

-the effective annihilation section of the p's in matter in 
the telescope. This ef fee tive ' section is parametrized by: 
a=48A 1 / 3 B - ’ 868 A=target mass and 8=v/c . In /figure 7 we see a 
comparison between this formula and the experimental results, 
-multiplicity of the annihilation pions (see figure 8). 

-the energy distribution of the annihilation -pions (figure 3) 
-the instaneous disintegration of the tt° into 2y . 

-the interaction of the y in the calorimeter producing a 
Compton e or an e+,e~ pair. 

The annihilation pions are isotropically generated in the 
center of mass, their direction and their energy being restored 
in the laboratory system by means of the Lorentz transform. The 
distribution of the pions in time is generated by the GENEV sub- 
program borrowed from the CERN-LIBRARY W-505 . For illustration. 
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we show in figure 9 the distribution of the effective masses 
7i+ tt + and n~ it - in the p-p annihilation at rest. The secondary 
particles emitted at the time of the annihilation are followed 
in the various detectors with integration of the Cerenkov light. 

At the time of an annihilation- in flight, the number of pions 
emitted is appreciably the same as for an annihilation at rest: 
only their kinetic energy increases. Figure 10 shows the distri- 
bution of the antiproton signals in the calorimeter, on the pes- 
simistic assumption that the calorimeter efficiency is 15 p.e/cm. 
The distribution width reflects the fact that certain secondary 
annihilation particles may escape from the calorimeter. We there- 
fore see that the annihilation of a p typically produces several 
thousands of photoelectrons. By comparison, a proton of velocity 
v/c=l passing through the calorimeter would produce about 1500 p.e 
(it would of course discard p as a candidate by the high signal 
it would produce in Cl) . The upper energy limit of the telescope 
is about 600 Mev/n which corresponds to the energy of a p (or a 
p) which is totally decelerated in the calorimeter. The maximum 
light emitted by a particle in the calorimeter would be achieved 
for a Deuterium of 600 Mev/n which would produce about 60 p.e in 
C2, or 25 times less than the mean antiproton signal. The anti- 
helium signal will be 4 times the p signal, i.e. an average of 
about 6000 p.e. 

The antimatter detection efficiency is determined by 2 factors 

-the proportion of antiparticles which become annihilated in 
module 1; 

-the proportion of antiparticles which become annihilated in 
the calorimeter, but one or several of whose secondary products 
rise in module 1. 

In the 2 cases, the signals of the various detectors of module 
1 will very probably be inconsistent and the event will be discard- 
ed. In figure 11 we showed the telescope efficiency of the p's 
calculated by the -simulation program . The efficiency is' lower 
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for i) the probability of annihilation in module 1 is greater and 
ii) the annihilation point in the calorimeter is closer to module /7 
1 and therefore the probability of having a secondary ascendant is 
greater. The 2cm thick lead plate is used to absorb part of the 
ascendants, the electrons in particular. 

In conclusion, the efficiency between 200 and 480 Mev/n is 28%, 
it is therefore 53% between 480 and 600 Mev/n. 

V - Background Noise 


We may consider various causes capable of simulating an anti- 
proton event (the problem is much less critical in the case of He). 

1 - The nuclear interactions in the calorimeter: the inter- 
action products are either slow nuclei (whose energy is quite a 
bit lower than the Cerenkov threshold of water), or pions. The 
production threshold of relativistic pions in a p-p or p-n reac- 
tion is about 600 Mev. The energy of the particles in the calor- 
imeter being less than 600 Mev/n, the nuclear interaction pions 
cannot give a large signal in C2. 

2 - The diffusion of particle s in a calorimeter PM. In this 
case only one PM would be affected, in contrast to the p signal 
where several PMs are affected. The Cerenkov radiators of the 
calorimeter are analyzed by 2 distinct channels corresponding to 
2 half-rings (1 PM in 2). We may thus discard the events where a 
single half -ring produces a signal. 

3 - The p's of atmospheric origin: 5g/cm2 of residual atmos- 

phere, the p flux produced by the interaction of the cosmic ra- 

diation with the atmosphere is p/p<<10 3 * * * 7 at low energy. This is 

due to the fact that the probability of producing an energy p < 

1 Gev is extremely small. 
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4 - Random coincidences : let us assume that during the 
the acquisition system analyzes a proton event, the 
relativistic particle passes through the calorimeter without 
passing through the coincidence. This particle would produce 
a fairly large quantity of light in C2 to simulate an antiproton 
signal. We may estimate the flux of relativistic particles pas- 
sing through C2 to be (J) = 5xl0 3 p/sec. For a typical dead time of 
x=3ysec the probability of random coincidence is: 

Pc— (<f)T) 2 /2 = 10 _4 ' /8 

Our objective being to achieve a sensitivity of 10 4 , we see 
that the rate of random coincidences might contaminate the meas- 
urement of the p/p ratio (although not all random coincidences 
can simulate an antiproton event) . 

The CRAC experiment will be equipped with a system for dis- 
carding random coincidences. We will use a rapid monitoring chan- 
nel designed to measure the time separating the arrival of two 
successive calorimeter pulses (analogous to an anti -pile up system 
used on an accelerator) . 

VI - Telescope Performances 
VI. I - Antimatter 

Geometric Factor = 420 cm2.sr 

Half-angle of Opening of telescope = 12 degrees 


Energy 

Proton 

Telescope 

a p p 

(Mev/n) 

Flux 

(p/sec) 

Efficiency 

Sensitivity 

(/hour) 

200-480 

24 

0.28 

4xl0~ 5 

480-600 

8 

0.53 

6xl0~ 5 
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We thus see that with a balloon flight of about thirty 
per hour, we achieve a sensitivity of about 8x10 7 in the de- 
termination of the p/p ratio. The sensitivity at the He is not 
as satisfactory, because the He interactions in module 1 will be 
greater and the mean number of secondary ascendants will be 4 times 
higher than for the p's. The telescope efficiency at the He's is 
5% and 15% respectively for the low and high energy ranges. This 
corresponds to a sensitivity of about 1X10 5 for the measurement 
of the He/He ratio. Conversely, we may consider expanding the 
energy range to 800-900 Mev/n for the detection of the He's, which 
would then give a sensitivity of about 5x10 4 , or one order of mag- 
nitude below the current upper limit. The situation for the p/p 
ratio is summarized in figure 12. One should note here that the 
measured energy range does not correspond to the energy of the 
cosmic radiuses in the interstellar environment due to the solar /9 
modulation. In the table below, we give an estimate of the number 
of detectable p's for 2 solar modulation conditions and a flight 
of 30 hours. 


presumed model 


Solar 
Activi ty 

Energy loss in 
the :solar cavity 
(Mev) 

primary p 
ref. (3) 

closed 

galaxy 

leaky 

box 

Minimum 

300 

210 p 

11 p 

2 p 

Max imum 

600 

4 0 p 

6 p 

1 P 


VI. II - Isotopes 

The CRAC telescope can separate the isotopes from the cosmic 

radiation in an energy range from 600 to 1000 Mev/n. The detec- 
tion principle is the following: a particle which triggers the 
Cerenkov Cl penetrates into the calorimeter C2 where it undergoes 
an energy loss per nucleon which is dependent upon the ratio Z 2 /A 
(Z-charge, A=mass) . In the calorimeter, the particle will emit 
Cerenkov light as long as its energy level is greater than 480 Mev/n. 


10 



Therefore, for an equal signal in Cl, a heavy particle will give 
a greater C2 signal than that of a lighter particle. In figure 
13 we see the C2 signal as a function of the Cl signal for hydro- 
gen aid deuterium. Figure 14 shows the mass spectra for the pre- 
ceding figure: the mass resolution is on the order of 0.12 mass 
unity. In this estimate, no correction was made of the angle of 
incidence (O<0<12 degrees) which would improve the resolution even 
more. 

VII - Specifications 


The CRAC experiment is planned to be launched by balloon. 

The telescope will be inside a pressurized pod. The problem of 
heat regulation is in large part solved due to the large quantity 
of water involved in the experiment. The experiment requires the 
use of 72 PM of photocathodes 120 mm in diameter for the Cerenkov 
and 80 PM 37 mm in diameter f or the hodoscope . The flight time 
will be equipped with 8 high-speed PM. 

Each Cerenkov and each scintillator layer has its own high- /10 
voltage power supply, i.e. 10 HV power supplies plus one HV for 
the flight time. 

Each hodoscope PM has its own analysis channel. Conversely, 
each Cerenkov only has 2 analysis channels: the 12 PM are assemb- 
led in groups of 6 by adding the signals of 1 PM to 2 . The Cer- 
enkov will be equipped with square root amplifers (as on HEAO) 
to accommodate enough dynamics for the isotopic analysis. For 
the hodoscope signals, one linear amplification is adequate, since 
the dynamic is only about 100 from Z = 1 to Z=8 . 

A microprocessor on board will have the function of i) sort- 
ing the events in real time and ii) adjusting, if necessary, cer- 
tain test parameters such as thresholds or acceptance criteria. 

A diagram of the collection system considered is shown in figure 
15. The Ranging rates is estimated at 20-30 kbits/sec. We are 
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also planning to install one (or several) magnetic recorders 
on board which would serve as a redundancy to the telemetry 
and which would make it possible to continue the acquisition 
in the event that the experiment is out of the receiving sta- 
tion's range. Such recorders now exist on the market with a 
capacity of 65 Moctets for a weight of 2kg. 

Telescope Weight 


Water 

Lead 

Scintillators 

PM 

Mechanical 

Electronics 

Batteries 

Pod 

Total 


550 kg 
115 kg 
4 0 kg 
12 0 kg 


1500 kg?? 


VIII - Implementing the Experiment 

The experiment shall be performed in 2 years in the labora- 
tory of the Saclay Astrophysics Service. A calibration of the 
experiment is planned at an altitude of about 3000 m. The Part- 
icles flux at this altitude, although lower than that found in 
the stratosphere, makes it possible to test a certain number of 
telescope functions such as anticoincidence and the anti -pile-up 
system. Furthermore, due to the fact that the antiproton flux 
is estimated at p/p=10 4 at an altitude of 3000m, we can check 
to what extent the background noise can simulate an antiproton 
annihilation. Finally, the D/H ratio being on the order of 10 3 , 
we may also measure the isotopic resolution of the experiment in 
conditions very similar to those of a stratospheric flight. 

The telescope will see about 3xl0 4 protons and 30 deutons 
per day between 400 and 1200 Mev/n. A calibration of some 15 
days at an altitude of 3000m should make it possible 'bo become 
sufficiently familiar with the instrument for a delicate interpre- 
tation of the future flight. The calibration data will aslo serve 


/II 
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to develop tifte'analysi s programs of the experiment. It would 
also be desirable to test the antiproton experiment with CERN 
to verify the parameters established by simulation, efficiency, 
p signal, etc . ) . 

The launching site under consideration is at Fort-Chunchi 11 
(Canada) where the geomagnetic cut-off is about 200 MV. The 
optimal period for the launching would be 1986-1987, coinciding 
with the next lowest solar activity (see figure 16). 

The p interaction distance in the atmosphere being about 
50 g/cm2, it is advantageous to fly at a residual pressure of 
a few mbars. To achieve this objective, we will use a 750000 m3 
balloon. 
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IX - Implementation Schedule and Financing 


/12 


1983 

1984 

1985 

1986 

BUDGET 

1983 

1984 

1985 

1986 
Total 


Physical, electronics, mechanical studies 
Flight Time Completion 
Cerenkov 1 Completion 

Hodoscope and anticoincidence mock ups 


Hodoscope completion 
Calorimeter completion 

Beginning of integration and ground testing 

Anticoincidence completion 
End of integration and ground testing 
Calibration at 3000m altitude 
p calibration at CERN (if feasible) 

Balloon flights in spring and/or autumn in Canada 


in thousands of TTC francs 


Flight time 

+Cerenkov 1 150 

Data acquisition system 150 

300 

5 Cerenkov 400 

Mechanics 150 

Electronics + Power 

supply 500 

Hodoscope 150 

1200 

Mechanics 180 

Anticoincidence 120 

Missions 50 

Recorders _50 

400 


Miss ion s+ transport 


(Canada) 250 

Batteries _50 

300 

Cost - - 2200 
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RATIO P/P 





Annihilation 
duction of c 
: S TT + IT - e + e 






1 ! I 






C 2-1 
C2.2 
C 2.3 
C2.4 
C2.5 



Figure 4 - Diagram of the CRAC telescope 
Module_l 

Tl TOF+co incidence scintillator 5 0 cm h=3 cm 

51 10 scintillator rods X view 50X5X3 cm 

52 ditto,. SI, y view , , _ 

Cl Cerenkov with water <j> = 60 cm h=15 cm 

53 ditto Si X view 

5 4 ditto SI Y view 

T2 ditto Tl 

Calorimeter 

C2-1 Cerenkov with water $=80 cm h=20 cm 
C2-2 ditto C2-1 

C2-3 ditto C 2-1 

C2-4 ditto C 2-1 

C 2-5 ditto C2-1 
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* C®**cnWov pulte httgM «rv farm* of t© 



P c 


Fig 5. Response of a Cerenkov 
with water as a function of the 
particle velocity. The curve 
represents the expected response 
without scintillation. The re- 
sidual light under nominal Ceren- 
kov is due to the 6 rays (title 
of reference (7)). 




S=S(pm)/ S(tot) 


Fig 6. Number of photoelectrons collected per 
cm of water for a particle Z=l,v/c=l as a func 
tion of the ratio S of photocathode surface to 
the total internal surface of Cerenkov. The 
experimental point is (taken from reference (6) 
The ratio S for the CRAC Cerenkov are also 
shown . 













1000 


2000 


3000 


N Photoelectrons 

Fig 10. Distribution of the light_collected 
in the calorimeter at the time of p annihil- 
ation with n=15p.e/cm. 
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Fig 11. Efficiency of the CRAC telescope 
for the detection of the p's as a function 
of energy. The full curve represents the 
fraction of incident p's which become an- 
nihilated in the calorimeter^ The dotted 
curve gives the fraction of p's annihilated 
in the calorimeter without secondary parti- 
cles rising in the telescope module. 

-a -Cerenkov threshold. 
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SIGNAL CERENKOV Cl 

Fig 13-Distribution of the Cl signals versus C2-3 
for hydrogen and deuterium between 600 and 700 Mev/n 
The calculation was performed with 50 p.e/cm for 2 
Cerenkovs. 
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HYDROGEN 


Fig- 14 -Mass spectrum of H and D obtained by 
adding the events according to the parallels 
shown in figure 13. The mass resolution is 
on the order of 0.12 mass units. 
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Key: Fig 15 -Diagram of the data acquisition system for the CRAC 

exper iment. 


1-Coincidence system 4x12 bits; 2-Flight time 1x12 bits.; 
3-Hodoscope 8x10 bits, 8x12 bits; 4 -Cerenkov channels; 

5 -Anti-stack-up system 4x8 bits; 6 -Management; 7-Format device; 
8 -Un switching; 9-//<;Ghartnel;' 10-Series channel; 11-Ground 
station HP 9816; 12-BUS Discris threshold 

Discris addresses 

Return of threshold measurements. 
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